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Abstract 

We show that the weight functional for color sources in the classical theory 
of the Color Glass Condensate includes a term which generates Odderon 
excitations. Remarkably, the classical origin of these excitations can be 
traced to the random walk of partons in the two dimensional space spanned 
by the SU(3) Casimirs. This term is naturally suppressed for a large 

nucleus at high energies. 
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1 Introduction 


Hadronic cross-sections have been successfully described for some time [I] 
in terms of non-perturbative t-channel exchanges known as Pomerons and 
Reggeons. In perturbative QCD (pQCD), the Pomeron arises from the color 
singlet exchange of two reggeized gluons in the t-channel [2]. This BFKL 
Pomeron, leads to cross-sections that grow rapidly as a power of the energy. 
Next-to-leading order corrections to the BFKL Pomeron |3j, as well as uni- 
tarization and saturation corrections, are topics of active interest in high 
energy QCD |4j. 

The Pomeron, which corresponds to a C-even exchange in the t-channel, 
has a C-odd partner dubbed the Odderon In pQCD, it arises from the 
color singlet exchange of three reggeized gluons in the t-channel [Ej. This 
BKP Odderon, in contrast to the BFKL Pomeron, has solutions that grow 
at most logarithmically as a function of the energy [7j. A nice review of the 
Odderon can be found in Ref. [HJ- 

Recently, Kovchegov et al. (KSW) |J . working in the dipole picture m 
EH, computed non-linear saturation corrections to the BKP Odderon equa¬ 
tion in the limit of large number of colors N c and large atomic number A. 
These are analogous to the saturation correction derived previously t .f2] for 
the BFKL Pomeron. These saturation effects modify the BKP result, leading 
to an Odderon contribution to high energy amplitudes that decreases with 
energy. 

Understanding the origins and role of the Odderon and its Pomeron 
partner, is essential for a consistent formulation of high energy QCD. One 
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such attempt is the Reggeon effective field theory of Lipatov na An al¬ 
ternative approach is the Color Glass Condensate f.Qj (CGC), where the 
degrees of freedom are parton fields at small x’s of interest coupled to static 
parton sources at larger values of x E3; their dynamical evolution with x is 
described by the JIMWLK Wilsonian renormalization group equations L®- 
Understanding the correspondence between the two approaches will help clar¬ 
ify our understanding of QCD at high energies. 

With this aim in mind, Hatta et al. HZ! identified the dipole and baryon 
operators with Odderon quantum numbers in the CGC framework. The 
JIMWLK equations for the evolution of these operators are shown to be 
equivalent to the BKP equation for low parton densities. In the large N c and 
large A mean held approximation, the KSW evolution equation for the Odd¬ 
eron dipole operator is recovered. The JIMWLK equations, in full generality, 
reveal a complex pattern of Pomeron-Odderon couplings under quantum evo¬ 
lution. Ref. HZ! establishes a specific correspondence between the framework 
of Reggeon held theory and the Color Glass Condensate. 

The authors of Ref. HZ! also argue that Odderon excitations are absent 
in the McLerran-Venugopalan (MV) model [^j , the classical effective theory 
(without quantum evolution) of the CGC. Their reasoning is correct for the 
original MV ansatz which assumed a Gaussian random distribution of the 
large x color sources in the CGC. However, we showed in a recent paper 
that the MV model for an SU(N C ) gauge theory where N c > 3 includes 
additional contributions to the distribution of sources ra. In this note, we 
will show that these contributions, sub-leading in the large A limit, generate 
the Odderon. For SU(3), this term in the path integral of the MV model is 
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proportional to the cubic Casimir operator. As in the small x limit, where 
(quantum) saturation suppresses Odderon excitations, (classical) saturation 
will be shown explicitly to suppress these as well. 

This letter is organized as follows. We begin by summarizing briefly 
the MV model emphasizing the coarse graining assumption implicit in the 
model. We state recursion relations (derived previously in Ref. |18j ) which 
describe the random walk in color space of SU(3) quark color sources. We 
extend the discussion of Ref. Pi to explicitly write down the path integral 
for the SU(3) MV model. With this path integral, we compute the Odderon 
operators identified in Ref. UH both in the weak field limit and to all orders 
in the classical field. In the weak field limit, for A = 1, our results for these 
agree with Ref. im 

2 The MV path integral for a large nucleus 
and an SU(3) random walk 

The MV model is a simple model which captures several key features of 
QCD at high parton densities P). The model can be generalized and the 
JIMWLK evolution equations we referred to previously arise from a proper 
accounting of quantum corrections in the Leading Logaritmic Approximation 
in x .ED.- The degrees of freedom are valence partons, which are random, 
static, light cone sources, and wee parton fields which couple coherently to 
valence sources all along the Lorentz contracted width of the nucleus. The 
kinematic condition for this coherence is x -C A -1 / 3 . How many random 
sources the wee parton fields couple to depends on the typical transverse 
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momentum of the wee parton 1 . A wee parton with momentum p± resolves 
in the transverse plane of the nucleus, an area (Aa;j_) 2 ~ 1 /p\. The number 
of valence partons it interacts simultaneously with is then 



( 1 ) 


which indeed is proportional to A 1 / 3 since A va i ence = 3 • A in QCD. This 
counting of color charges is only valid as long as p± > Aq CD . 

The classical effective Lagrangean for the MV-model, formulated in the 
infinite momentum frame (P + —> oo) and light cone gauge (A + = 0) has the 
form, 


£ = S[A,p]+iW\p]. 


( 2 ) 


Here S[A, p\ has two terms. The first is the usual QCD Field Strength 
tensor squared describing the dynamics of the wee parton gauge fields AC 


The second describes the coupling of these fields to the source color charge 


density p. These have been discussed previously |2j and will not concern us 
further here. 

Our focus will be on the final term-corresponding to a Boltzmann-like 
weight exp {—W[p]) in the generating functional for the effective action. It 
denotes the likelihood of different p configurations. In the original MV model, 
the weight functional for a very large nucleus had the Gaussian form 



(3) 


lr The wee parton is soft only in its longitudinal momentum-its transverse momentum 
may be large. 
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where the average color charge squared per unit area per color degree of 
freedom, n\ is simply determined to be mm, 


2 


(4) 


It is of order A 1 ^ fm -2 . For a large nucleus, this scale is the only dimensionful 
scale in Eq. El Thus chs(ha) ^ 1, an d all features of the theory can be 
computed in weak coupling. 

In Ref. [IB] , it was shown explicitly that the assumption about the Gaus¬ 
sian distribution of charges for k » 1 is exact only for quark sources in an 
SU(2) gauge theory 2 . For SU(3) quarks, there is an additional contribution 
to Eq. 0 Since the argument has been stated at length previously in Ref. (IS], 
we shall only outline key features below with emphasis on novel results. 

The problem of interest can be formulated as follows. Given k random 
577(3) quarks in a coarse grained cell of transverse area (Ax_l) 2 , what is 
the distribution of color representations? This problem can be formulated 
straightforwardly as a recursion problem. 

When one adds a quark to an (m, n) representation of SU(3), the result 
is a sum of three irreducible representations 3 , 


(1, 0) X (ra, n ) = (m + 1, n) + (m — 1, n + 1) + (ra, n — 1) (5) 

2 It is also true for gluon and quark-anti-quark sources in an SU(3) gauge theory. For 
a detailed discussion, see Ref. ms. This fact also explains why only the Gaussian term 
was obtained in the computation of Ref. EDI, since there the nucleon was modeled as a 
quark anti-quark dipole. 

Representations in SU( 3) can be expressed as (m,n), where m and n are integers 
which label the upper and lower tensor indices respectively. A fundamental quark 3 state 
in 517(3) is (1, 0) while an anti-quark in the 3 state is (0,1). 
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One can deduce therefore that the multiplicity of a particular represen¬ 
tation iV^ 1 ) in the distribution of k + 1 quarks (among all possible repre¬ 
sentations) is given by the recursion relation 


l\fi k + 1) — ]\f ( k ) /\d fc ) i /\r( fc ) 

(k) (k) 

for m, n > 0 with the boundary condition N^-i = AT/ n = 0 and the initial 
condition, Aq'q = 1 e/se -/V® n = 0 . 

As shown in Ref. [THj, the multiplicity Nff n can be expressed in terms 
of the tri-nomial function Gk- m ,n defined as 


G k: 


m,n 


One hnds that 


k\ 


( k+2m+n \ | / k—m+n \ | ( k—m—2n \ | 
V 3 )■{ 3 M 3 )■ 


(7) 


7 Ti,n m,n Gk: m+3,n “1“ Gk: m,n+3 

Gk : m+2,n—1 G^/c: m— l,n+2 Gfc m+2,n+2 (8) 

satisfies the recursion relation Eq. © with the stated initial and boundary 
conditions. 

The two Casimirs for a given SU(3) m,n state are the quadratic and 
cubic Casimirs [23], 

Do 71 ’"' 1 = — ( m 2 + mn + n 2 ) + (m + n) 

_j_ 2 n + 3 ) ( n + 2m + 3) (m — n) . (9) 
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Note that while the former is symmetric in m and n, the latter is anti¬ 
symmetric under exchange of m and n. For large k, one can use the Stirling 
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formula to simplify Eq. Q and further use Eq. El to obtain the very simple 
result: 

G k :m,n ~ 777 — exp (-3D™’ n /k) (l + 3D™’ n /k 2 ) . (10) 

ZK'K v ' 

The multiplicity, defined in terms of the G’s in Eq. |S1 can be approxi¬ 
mated in the large k limit as 


« 2 d 3 m G k: m,n + 2d 3 n G k:m , n -3 d m d 2 n G k: m , n -3 d n d 2 m G k:m , r 
27 m n (m + n ) 37 +fc 


k 3 


2 A) 7T 


exp (-3 D?’ n /k) (l + 3 D™' n /k 2 ) . (11) 


While Eq. HH1 was derived in Ref. pS], the final expression in Eq. is new. 
The fact that the correction due to the cubic Casimir in Eq. is identical 
to that in Eq. [TJHis non-trivial and remarkable. From Eq. El one can deduce 
that the peak of the distribution is at m = 0{y/k) and n = 0(Vk). Cor¬ 
respondingly, the cubic Casimir introduces a correction of size 0(1/ \/k) for 
large k. 

The dimension of an SU(3) representation is 

(m + l)(n + l)(m + n + 2) mn(m + n) 


dmn 


( 12 ) 


The probability to hnd the system of k quarks in an (m, n) representation 
is therefore proportional to d mn N^\. One can then define a normalized 
distribution 1 in the sense that 

1 = [ dmdnm 2 n 2 (rn + n) e~ 3L>2 ^ k (l + 3D™’ n /k 2 ) . (13) 

_ 4fu 7T Jo _ v 7 

4 One should not interpret the integrand as a probability density when m and n are 
large compared to \/k. But this is academic since the probability to find such large to’s 
and n 's is exponentially suppressed. 









Since the cubic Casimir changes sign if m and n are exchanged, the correction 
term gives a vanishing contribution to the integral. 

In Ref. P3j, we proved that the right hand side of Eq. El] can be re¬ 
written exactly as an integral over classical color charges 5 * , 

1 ~ J d 8 Q e -^Q 2 A+3^3(Q)/A: 2 (14 ) 

where the cubic Casimir term in the integral is interpreted as a perturba¬ 
tion to the dominant quadratic Casimir term. Here Q = (Qi, Q 2 , ■ • •, Qs) 
is a classical color charge vector defined by |Q| = RQ a Q a = D™’ n and 
Qii" m 1 Q& are its eight components. The cubic Casimir, in terms of the 
color charges, is D 3 (Q) = dabcQ a Q b Q c , where d a b c is the symmetric tensor in 
the SU( 3) Lie algebra. 

Q a is the classical color charge in a single coarse grained box. One can 
define a local classical color charge density p a in terms of Q a : 


Qa = 


f d 2 xp a (x) « i(Ax_L) 2 p a (x). 
J(Ax 1 I 2 Q 


(15) 


9 J( l*i)’ 9 

With this definition, and k defined as in Eq. E one can write down a path 
integral for the distribution of the color charges as 

Pa (x)Pa (x) d a6c p a (x)/(x)p C (x) 


exp 


gLx 


2/4 


k a 


(16) 


where /i 2 4 = g 2 A/2nR 2 A was defined in Eq. E For the weight of the cubic 
Casimir, we obtain, 

g 3 A 2 N c 


k a = 


7 t 2 R 4 


(17) 


As we will now see, the cubic term in Eq. EI]is precisely the term that gives 
rise to the Odderon. 

5 Our derivation did not include the cubic Casimir but that does not change anything 

since that term integrates to zero. 
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3 Odderon in the classical effective theory 

In Ref. [T7], the authors identified the Odderon exchange operators in two 
situations, a) a dipole scattering off the CGC, and b) a three quark baryon 
system scattering off the CGC. We first consider the former where the C-odd 
“dipole odderon operator” is defined to be 

0(x, y) = ^Tr(vJV„-VjV I ) , (18) 

where 

Vj = ^-oo,oo(x) = P exp (ig j dx~a a (x~,x.)t a ^j , (19) 

and VIqq oo = Voo-oo. Here t a is a generator in the fundamental representa¬ 
tion of SU( 3) and a a is the CGC gauge field. In the weak field limit, one 
can expand out the path ordered exponentials-the lowest order contributing 
to the Odderon is of cubic order in a Eh 

0(X, y) « - oQ (“1 - <) K - . (20) 

Now, from the classical equations of motion for gauge fields in the CGC, the 
gauge field can be expressed in term of the classical color charge density p a 
in Lorentz gauge as 

a * = hf d2z ln ( (x-W ) pa(z) ’ (21) 

where A is an infra-red cut-off which will not appear in the final result. 
From this relation therefore, one expects O oc d abc p a p b p c . In the original MV 
ansatz, the expectation value of this Odderon operator would be zero. 
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However, from Eq. 1161 we now know that there is indeed a contribution 
(albeit suppressed) in the source term of the MV model which has the color 
structure of the cubic Casimir operator d abc p a p b p c . One can therefore com¬ 
pute the initial condition for quantum evolution of the CGC directly from 
the MV model. 

To evaluate ( O ), with the measure in Eq. [THl one needs to compute the 
integral 

X = dabcd ~ alz f[dp]e~f d2xp2/2 ^ 

Zq ka J 

x j ifyf>°(y)p\y)p c (y)p‘(u)p\-v)?(v,). ( 22 ) 

As prescribed previously, we have expanded out the cubic piece in Eq. fTTH 
here. Z 0 corresponds to this measure. 1 can be easily computed using 
Wick’s theorem. Ignoring disconnected pieces for the moment, one obtains, 

1 = 6^4m r ,# 2 y (p a (y)p s (u))(p b (y )p^(v)) (p c (y )p e (w)) (23) 

Zo ka J 

where we used the fact that d abc is a totally symmetric tensor. In the MV 
model, the 2-point correlator is simply 

(p“(x)A y))=^*< a) (x-y). (24) 

Using Eos. UHlIiHl to compute the expectation value of O, one obtains, 

(C>(x, y)) = 9 V A d abc d abc f d 2 u (G(x - u) - G(y - u )) 3 

4 N c k a J 


= a 


I M-iiWHUvs r d > u 

4 7T Tn IVj? J 


In 3 


x — u 


u 


(25) 
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where we defined G(x — y) = In ( ^ x _ ' 2A2 ) /4vr. Also, note that and ka 
were defined in Eq. E]and Eq. Irrespectively, ro = 1.12 fm and we used the 
identity d abc d abc = (TV 2 - 1) (TV 2 - 4)/TV c . 

In performing this computation, we ignored the disconnected pieces 
coming from the Wick expansion of Eq. These pieces vanish because 
they involve a summation with two contracted indices of d a bc- 

We now turn to the Odderon operator corresponding to a three quark 
baryon state scattering off the CGC. In that case, the relevant operator 

is pill ITT ] 

B(x - y - z) = 3k V (x) V»(y) ~ h - c ) ■ (») 

where VJ is defined in Eq. E3 Again, expanding out to lowest non-trivial 
order in a, the expression for B(x,y,z)in the weak field limit is 

« ^ <1* {(«; + a; - 2a3 (aj + 2^) « + ^ - 2«;)) . (27) 


Replacing a —> p using Eq. EH using Wick’s theorem and computing the 
average over color sources using Eq. EH one obtains, 

6 


(6) = 


d u (G(x — u) + G(y — u) — 2 G( z — u)) 


x (G(z — u) + G(x — u) — 2 G(y — u)) 
x (G(y — u) + G(z — u) — 2 G(x — u)) 

[*uln 

7T Tn TVr 2 

x In 


ai 
24 7T) 


x - 

- u 

|y 

— u 


z — 

■ u 

2 


|y - 

- u 

z 

— u 

1 In 

z - 

- u 

X 

— u 


x — 

■ u 

2 

|y - 

- u 

2 


(28) 

y - u r 

For TV C = 3 and A = 1, both Eq. EH and Eq. EH can be shown to agree 
exactly with Eqs. 5.5 and 6.2 respectively in Ref. El- To compare, we need 
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to assume the nucleon is dilute and the fact that the result of Ref. DU is for 
a single quark (namely divide our result by 3). In Eq. Eland Eq. El the 
factor of A 1 / 3 accounts for the possibility that the Odderon couples to A l A 
nucleons in the nucleus. 


4 Result to all orders in the gluon density 


We shall now explicitly evaluate the expectation value of the dipole odderon 
operator in Eq. El to all orders in g 2 p\/k A . In performing this computation, 
we have to take into account the path ordering of the gauge held in the x~ 
direction-see Eq. El The expectation value of 0(x, y) is given by 

(O) = J du~ J d 2 u^p a (u~,u)p b (u~,u)p b (u~,u) 

x Tr (V(x)V(y) - E(y) f E(x)) ) 

where the average is now done with the “smeared” Gaussian weight 

Pa(z-,z)p a (z-,z)\ 


W[p\ = exp i — Jdz J d 2 z 


(29) 


(30) 


2 d 2 (z~) J 

Note that (p a (x~, x)p b (y~, y)^> = p 2 (x~) S ab 5(x~ — y~) 5^(x — y). To eval¬ 
uate (O), we use the following identity, 


d 


6 -W\P\ = - Pa{u)p b (u)p c {u) W[ P \ (31) 


abc 5p a (u) 5p b (u) ~5p c (u) " Lrj p 6 (u )' 
using the fact that IE is a Gaussian and d aa b, summed over a is zero. Re¬ 
placing the LHS of this identity, which appears in Eq. El with the RHS and 
integrating by parts, we can write O as the sum of two terms, where the first 

is 

dabc f i — r / f I & & & 


(O), = 


j du p 6 (u ) J d 2 u^- 


2 iN c K A J ~ r ~ v '~ 'J ™ x dp a (u) 5p b (u) 5p c (u) 

x Tr (E_oo,oo(y)Ko -oo(x)) ) • (32) 
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The second term (0) 2 is obtained by changing the overall sign in (0) 1 and 
exchanging x and y. After some algebra 6 , 


(O) i = 


£ 3 (N 2 - 4) 
2iNhiA 


J du /i 6 (w ) J d 2 u ^ c 4iV (G 3 (x - u) - G 3 (y - u)) 


x 


J 1 00 ’ 00 + y (<T 2 (y - u)G(x - u) - G(y - u)G 2 (x - u)) • X 2 °°’~ 


(33) 


where X“’ 00 = -oo(x)) and 

XT’" 00 = (4K-,oo(y)f4o, u -(x)4K-,-oo(x)lA 00i „-(y)^. X 2 can be simplified 
using the locality of the distributions in x~ . Further, since 


X V ,U ( 

x ’ = exp ( — 


9 2 t a t a 


dz fi 2 {z ) / crz (G(x — z) — G{y — z))' 


and t a t a oc I 3 X 3 , it commutes with all t a , and the two sets of (V^V) in the 
second term of Eq. EH] can be combined to form Substituting these 

results in Eq. 021 and combining it with (Oo) (obtained by exchanging x and 
y), we have, 


(o) = 4 

x Tr exp 


(^- 4 )(^- l) ^ 1/3 


4 7 r Tq N* 

g 2 l-l 2 A tata 


d 2 u In 3 


u 


u 


(fz (G(x - z) - G(y - z))' 


(34) 


We have replaced here f^dz - g 2 (z~) = fi\. Note that Tr exp() — N c — 

{{N 2 — l)/4 - - -, so the first term in the expansion is identical to the weak 

6 We have used the identities d a bct a tb = \ Nc N ~ A t c and d a b c t a tbtc = — -j -1) 1 where 
1 is the 3x3 identity matrix. We have also used the cyclic property of the trace. 

' The local structure in x~ (see the expression after Eq. ESI ensures that only adjacent 
Wilson lines can contracted |25b This simplifies X“’ _oc to read 

X2 ta (Ez - jQO (y(A)} ta ( Ez — , — oo (x) ^—oo,u~ (y)} 


14 













field result in Eq. [251 As discussed previously, the argument of the expo¬ 
nential here can be expressed in terms of the saturation scale Q 2 A of a large 
nucleus EH [251 d. Our expression in Eo. 1521 agrees with the initial condition 
(modulo geometrical factors) for Odderon evolution suggested in Ref. || in 
a diagrammatic approach. In the large A limit, when the saturation scale 
is large, Eq. [52] tells us that the Odderon contribution will be suppressed. 
Of course, there is a further suppression that comes in from quantum evolu¬ 
tion [7 . 

We shall now outline the computation of the Baryon Odderon operator 
to all orders in the classical color held. The full computation, while perfectly 
feasible and straightforward, is extremely tedious for reasons which will be¬ 
come clear. This will be left for a later date. We begin with an alternative 
expression to that in Eq. [251 also given in Ref. na, which will suit us better: 

B(x,y,z) = W [Tr(V?Vy Tr(gVj - Tr^V^Vj - h.c.' 

= /Cl - /C 2 - h.c. (35) 

To compute the expectation value of this object, we follow the same pro¬ 
cedure as that for the dipole Odderon operator. Consider for example the 
expectation value of the Erst term. Using the identity in Eq. [^ we have, 

{K) x = [ du~ [ d 2 u( . d d d. . 

12 ik a J J \dp a (u)dp h {u)dp c [u) 

(V r _ 00i00 (x)U 00 _ 00 (z)) (U-oo,oo(y)Ko,-oo(z))). (36) 

One has an analogous expression for (/C 2 ). Unlike Eq. 02] where one had to 
manipulate 8 terms (2 from each differentiation), one now has 4 3 = 64 terms 
from (/Ci) and an equal number from (/C 2 ) to make a total of 128 terms! These 
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can be combined using Fierz identities eventually into Gaussian correlators 
of 2-point, 3-point and 4-point Wilson lines. A general technique to compute 
these was discussed in Appendix A of Ref. [26lj . 

5 Summary 

We showed in this note that, for SU(3) quark sources, the MV model in¬ 
cludes a term proportional to the cubic Casimir in the path integral. We 
obtained in Eq. El] a quantitative expression for this correction. Remarkably, 
this term derived from a random walk analysis of classical color charges is 
responsible for generating Odderon excitations in the classical effective the¬ 
ory of the Color Glass Condensate. Our results for the expectation values of 
the dipole Odderon operator and the Baryon Odderon operator in the CGC 
background field agree exactly with those of Ref. in the weak field and 
A = 1 limit. In addition, we computed the expectation value of the dipole 
Odderon operator to all orders in the classical field and showed that the result 
agreed with Ref. |5j. We outlined the computation of the Baryon Odderon 
operator-the detailed computation will be left for a future study. Our results 
suggest that Odderon excitations, while quantitatively suppressed, arise nat¬ 
urally in a consistent mathematical framework and are therefore relevant for 
a systematic study of high energy QCD. 
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